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The excitation function for the fusion-evaporation reaction ''^Ni + ^""Mo has been measured down 
to a cross-section of ~ 5 nb. Extensive coupled-channels calculations have been performed, which 
cannot reproduce the steep fall-off of the excitation function at extreme sub-barrier energies. Thus, 
this system exhibits a hindrance for fusion, a phenomenon that has been discovered only recently. 
In the S'-factor representation introduced to quantify the hindrance, a maximum is observed at 
i?s=120.6 MeV, which corresponds to 90% of the reference energy El"^^ , a value expected from 
systematics of closed-shell systems. A systematic analysis of Ni-induced fusion reactions leading to 
compound nuclei with mass A= 100-200 is presented in order to explore a possible dependence of 
the fusion hindrance on nuclear structure. 
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I. INTRODUCTION 

Heavy-ion induced fusion reactions have been studied 
extensively for more than forty years, especially since the 
discovery of the sub-barrier enhancement phenomenon 

IE ^ IS ' Coupled-channels descriptions were shown 
to explain the phenomenon successfully 0, • Recently, 
evidence was found for a strong hindrance of the fu- 
sion process at extreme sub-barrier energies, an effect 
for which there is no satisfactory explanation in present 
model calculations |1|. A systematic survey of existing 
data from the literature shows that there is a regular 
pattern to the energy dependence for the appearance of 
fusion hindrance For stiff, closed-shell colliding sys- 
tems, a significant maximum in the S'-factor is present 
as a function of the beam energy. This maximum signals 
the onset of sub-barrier hindrance, and the energy. Eg , at 
which it is located can be described well by an empirical 
formula 0: 

El'f = 0.356 {ZiZ2^)^ (MeV), (1) 

where ^ = AiA2/{Ai + A2). For softer systems, this 
formula provides an upper limit for the energy at which 
the 5'-factor has its maximum. 

Fusion hindrance at extreme sub-barrier energies is rel- 
evant, not only for understanding the dynamics of reac- 
tions between complex systems, but also for astrophysics 
and for the synthesis of superheavy elements. 
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The influence of nuclear structure on this hindrance 
behavior was first studied in a detailed comparison [loj | 
for the coUiding systems: ^^Ni + ^^Wi ^^^m + 
[13, 58Ni + 64Ni [13, and ^^Ni + ^^m^, where the 
systems are arranged in order of decreasing stiffness. For 
^^Ni -I- ^''Ni, an open-shell colliding system, the measured 
location in energy of the S'-factor maximum is about 9% 
lower than the value expected from Eq. ^ ^3 j while the 
other systems are more in line with the systematics. 

The aim of the present paper is to further investigate 
the hindrance phenomenon by measuring fusion evapora- 
tion for the open-shell system ^"^Ni -I- ^"^Mo. Compared 
to ^"'Ni, ^""Mo is a transitional nucleus with two protons 
outside the closed proton shell. Two earlier measure- 
ments of fusion excitation functions for ^^Ni -I- ^^^Mo 
can be found in the literature [isl Il6j |. The minimum 
cross section measured in these two experiments is about 
0.4 mb. Since the hindrance behavior is expected (from 
Eq. ^ to occur at much lower energies, the main aim of 
the present measurement was an extension of the excita- 
tion function into the nb region in order to localize and 
quantify this hindrance. 



II. EXPERIMENTAL PROCEDURE AND 
RESULTS 

The experiment was performed with ^^Ni beams in the 
energy range of 196-262 MeV from the superconducting 
linear accelerator ATLAS at Argonne National Labora- 
tory. The maximum beam current used was ~ 60 pnA. 
The high melting point of the target material, metal- 
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lie molybdenum evaporated on a 40/ig/cm^ carbon foil, 
prevented damage to the target by the relatively high 
beam current. The target thickness was constant during 
the experiment, as monitored with Si detectors. Thin 
targets with thicknesses of 8 or 18/j,g/cm^ were used 
in order to reduce the correction for target thickness 
in the energy regime of steep fall-off of the excitation 
fmiction. The isotopic abundance of ^"°Mo was 97.42%, 
with the remainder coming from ^*Mo (0.96%), ^^Mo 
(0.28%), 96Mo (0.34%), s^Mo (0.29%), ^^Mo (0.18%), 
and ^^Mo (0.53%). The selection of the beam and target 
combination excluded the possibility of background from 
fusion reactions coming from beam or target contami- 
nants Two surface-barrier Si detectors, located at 
± 43° with respect to the beam direction, served as mon- 
itors. The absolute cross sections for fusion-evaporation 
were determined by using elastic scattering measured 
with the monitors. 

The evaporation residues were identified and measured 
with the Fragment Mass Analyzer (FMA) 0| , which has 
been upgraded with the installation of a split-anode in 
the first electric-dipole 18] . The background, originating 
mostly from scattered beam, was greatly suppressed af- 
ter this upgrade. In order to push the cross section mea- 
surements to the lowest level, a new focal-plane detector 
with the configuration PGAC-TIC-PGAC-TIC-PGAC- 
IC, was used in the experiment. Here, the symbols PGAC 
stand for x-y position sensitive parallel grid avalanche 
counters |19| |. TIC for transmission ionization chambers 
|20j |. and IC for a large volume multi-anode ionization 
chamber. The first PGAC was mounted at the horizon- 
tal (x-direction) focal-plane of the FMA. 

The y-focus condition occurs about 70 cm downstream 
of the focal-plane, which is nearly at the middle of the 
last ionization chamber. The three sets of position sig- 
nals xi,yi,X2,y2,X3 and 2/3 were measured with the 
three PGAC's. The flight-times t2 and obtained from 
PGAC1-PGAC2 and PGAC1-PGAC3 were also recorded. 
Seven AE signals were measured with the ionization 
chambers (the first four AE1-AE4 in the two transmis- 
sion ionization chambers and the last three AE^-AE^ in 
the final ionization chamber, IC). Another time-of- flight 
signal, trf, measured the time difference between the rf- 
system of the accelerator and the first PGAC. Three ad- 
ditional AE signals from the PGAC's were also recorded. 
The three PGAC's and two TIC's operated at a fixed 
pressure of 3 torr of isobutane. The last ionization cham- 
ber had an adjustable pressure of 22 - 30 torr of isobu- 
tane. There was one 0.13 mg/cm^ Mylar pressure foil in 
front of the first PGAC, and another one with a thickness 
of 0.22 mg/cm^ located between the third PGAC and the 
last ionization chamber. This setup allowed for full track- 
ing of each particle detected in the system, and provided 
very good separation between evaporation residues and 
background events. A more detailed description of the 
new detector system will be published elsewhere pll |. 

For most settings of the FMA, two charge states of the 
residues were collected simultaneously. For the energies 
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FIG. 1: Two-dimensional plots of AEa vs. ts for the new 
focal plane detector system at the FMA, obtained at sub- 
barrier incident energies of 198.8 MeV (a) and 196.0 MeV 
(b). The isolated group in (a) (open circles) originates from 
evaporation residues, whereas the other events are caused by 
background. At E=196.0 MeV (b), 23 events (open circles) 
fall inside the acceptance window for fusion-evaporation. 



Slab = 260.5, 245.8, 209.1, 207.1 and 202.2 MeV, full 
charge state distributions were measured, while for most 
other energies, two FMA settings, i.e. four charge states 
were recorded. At the four lowest energies, only two 
charge states were measured. From the full charge state 
distributions, charge state fractions were determined for 
extrapolation to all other energies. The energy distribu- 
tions and angular distributions of evaporation residues 
were calculated with the statistical code PACE ,2^] . To- 
tal angular distributions for fusion-evaporation have been 
measured in Ref. for ^^Ni -I- ^"^Mo. In order to 

check the PACE code, calculations were compared with 
these experimental angular distributions of Ref. [Tsj ] and 
good agreement was found as long as the total calcu- 
lated angular distributions were taken as a weighted sum 
of the angular distributions of the different masses from 
our m/q measurements and folded with multiple scat- 
tering. Whereas in Ref. rather thick targets were 
used, the corrections from multiple scattering are small 
in the present experiment. The transport efficiencies of 
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the FMA were calculated with these angular distributions 
together with Monte Carlo simulations, using a modified 
version of the GIOS Code 23]. The large momentum 
acceptance, ± 10 %, and the large angular acceptance, 
6iab < 2.3°, of the FMA result in a high detection effi- 
ciency for the residues. 
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FIG. 2: Two-dimensional plots of AE5 vs. for the focal 
plane detector system at the FMA, obtained at sub-barrier 
incident energies of 198.8 MeV (a) and 196.0 MeV (b). The 
isolated group in (a) (open circles) originates from evapora- 
tion residues. The other events are caused by background. 
At 196.0 MeV, only one of the candidate evaporation residue 
events falls within the AE5 range. See text for details. 

The flight-time ts between PGACi and PGAC3, to- 
gether with a AE signal measured in a counter far behind 
the focal plane (e.g. AEq), was found to give the best 
separation of the residues from background events. Two- 
dimensional plots of AEq vs. ^3 are shown in Fig. ^ 
indicating the excellent separation achieved with these 
two quantities. Down to a cross section level of ~ 300 
nb, the evaporation residues could be identified on the 
basis of these two parameters alone. The events that fall 
within the expected window for evaporation residues are 
shown as open circles in Fig. QJl. Events shown by small 
points originate from beam particles, scattered from the 
beam pipe or from the area of the "beam stop" at the 
first split anode. The events in Fig. ^ correspond to 



a cross section of 242 nb obtained in a 12 hours run. A 
two-dimensional plot of AE<^ vs. ia, corresponding to the 
same events is shown in Fig.|2K- At this energy evapora- 
tion residues (open circles) are again well separated from 
background events (small points). 

At the lowest energy, i?beam=196.0 MeV, we find that 
23 events fall within the expected evaporation residue 
window as illustrated in Fig.^^- Most of these events are 
suspiciously close to an intense band of scattered beam 
particles which borders the fusion evaporation residue 
window. By examining the Ai?5 vs. ^3 spectrum shown 
in Fig.lSb, it is evident that only one of these events falls 
within the expected range of the AEc, signal. However, 
this event does not have the correct m/g-value as will be 
demonstrated below. 

An m/q spectrum containing the 23 events that fall 
inside the AEq — window is given in Fig. This 
spectrum is spread out over the whole m/q range and 
does not appear to be associated with fusion evaporation 
residue events seen at higher energies to be concentrated 
near channels 103 and 145 (see Fig. E^). The only can- 
didate event at i5=196.0 MeV, which satisfies the condi- 
tions for the AE^, Ai?g, and signals, appears in chan- 
nel 114 as seen in Fig. ^p. By further checking all the 
detector signals for this event, it was concluded that it 
is a background event, probably arising from scatterings 
of a beam particle at several locations through the spec- 
trometer. Hence, only an upper limit for the measured 
cross section could be determined at this beam energy. 
These results show that a background suppression fac- 
tor of about 4x 10^^^ can be achieved with the upgraded 
FMA and the present detector system. 

The uncertainties in the reported evaporation cross 
sections arise mainly from corrections for the charged 
state distribution, the FMA transport efficiency, the de- 
tector efficiency, and counting statistics. The total un- 
certainties for the evaporation cross sections are around 
10-17%, except for the measurements at the two lowest 
energies where only upper limits can be given. The cross 
sections are listed in Table L 

The fusion-fission cross sections for ^^Ni -I- ^"^Mo have 
not been measured previously. They were, however, 
calculated in Ref. d5:j with the code CASCADE pi ]. 
Similar calculations with the same parameters were per- 
formed for the present experiment. The total fusion 
cross sections are also listed in Table I, together with the 
fusion-fission cross sections. Rather large uncertainties 
were given for the fusion-fission contributions resulting 
in somewhat larger errors for the total cross sections at 
the highest beam energies. 

The experimental results for the total fusion cross sec- 
tions, spanning eight orders of magnitude, are presented 
as a function of laboratory energy in Fig. 01 (solid cir- 
cles). The incident energies have been corrected for tar- 
get thickness and for the steep energy dependence; these 
corrections are small because rather thin targets were 
used. For the lowest two energies no evaporation residue 
was observed. The results are shown as upper limits cor- 
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TABLE I: Measured fusion evaporation cross sections and 
calculated fusion-fission cross sections (using the code CAS- 
CADE) for the ''''Ni -I- ^°°Mo system. Large uncertainties 
were assigned to the calculated fusion-fission cross sections, 
leading to rather large errors in the total cross sections at the 
highest beam energies. A'^, is the number of charge states for 
evaporation residues which were measured in the experiment. 
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FIG. 3: Position {m/q) spectra obtained at energies of 198.8 
MeV (a) and 196.0 MeV (b). These events faU inside the 
windows in Fig. ^ or respectively. The event that also 
satisfies the AiSs condition is marked solid grey in panels b 
and c. See text for details. 



responding to one count in each case. The two earher 
measurements are also shown in Fig. 0] as open 

circles and open squares, respectively. Statistical errors 
are smaller than the size of the symbols. The three mea- 
surements are generally in good agreement within the 
quoted uncertainties, except for the lowest energy point 
ofRef. O. 



TABLE IL Structure input of low- lying states in Ni and 
^°°Mo. For ''"Ni, the B(EA)-values for the quadrupole tran- 
sitions are from |2rL l26|. the octupole strength is from Ref. 
[27|. The ^""Mo input is from Wi- 
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III. COUPLED-CHANNELS CALCULATIONS 



The coupled-channels calculations of Ref. repro- 
duced the excitation function quite well down to the 0.5 
mb level. It was found that with the same coupling pa- 
rameters, the calculations cannot reproduce the new data 
at lower beam energies. In fact, in the literature there 
are not many such calculations that can reproduce suc- 
cessfully experimental data for heavy systems. Since the 
nucleus -'^°'^Mo is rather soft, multi-phonon states and 
large coupling effects should be included in the calcu- 
lations and this represents a challenge. In the following, 
two sets of coupled-channels calculations are presented, 
including two- and three-phonon states, respectively. 

The nuclear structure input for the calculations is given 



TABLE III: Energies and reduced transition probabilities of 
the I = 0,2 and 4 transitions. 
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FIG. 4: Fusion excitation function for the system ®*Ni + 
^""Mo (solid circles). Statistical errors are smaller than the 
symbol size. For the two lowest energies, no residue events 
have been observed. The upper limits correspond to one 
count. Included in the figure are two previous measurements 
from Ref.|l5j| (open circles) and Ref. [l^ (open squares). 



in Table 11. The structure input for ^''Ni is the same 
as used in Ref. To] to analyze the ^''Ni + ^"'Ni fusion- 
evaporation excitation function. The nuclear quadrupole 
coupling, (3^ , for ^"^Ni was set ~ 10% higher than the 
value for Coulomb coupling, /3^, (for reasons discussed 
in Ref. [l^). For simplicity, it was assumed that the 
quadrupole mode corresponds to a perfect vibrator. This 
is not unreasonable, considering the uncertainty of the 
predicted one- and two-phonon couplings shown in Table 
111. The structure input for ^'^"Mo is also given in Table 
11. Here, the nuclear /3- values were also set ~ 10% higher 
than the Coulomb /3-values. 

The two-phonon (2ph) calculations include the 2+ and 
3~ one-phonon states in both nuclei, as well as mutual 
excitations and the two-phonon states listed in Table 11. 
The two-phonon octupole excitation of ^°°Mo was also 
included. The two-phonon states were treated as per- 
fect vibrational states. This gives a total number of 14 
coupled channels. The three-phonon (3ph) calculations 
additionally include all mutual excitations of the states 
considered in the 2ph calculations up to two mutual two- 
phonon states. Moreover, the estimated three-phonon 
excitation of the 2+ state in ^""Mo shown in table 11 was 
also included. However, states above 5.5 MeV were not 
included, so the total number of channels is 31. The ion- 



ion potential was parametrized as a Woods-Saxon well 
with a depth of Vo = -82.9 MeV, diffuseness a = 0.686 
fm, and nuclear radius Rn = 10.19 + Ai? fm. 
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FIG. 5; Fusion excitation function for "''Ni + ^°°Mo com- 
pared with several calculations described in the text. 

If the value of Ai?=0.06 fm was used in the two-phonon 
calculations, the result is similar to the calculations in 
Ref. 15], and it agrees quite well with the data for 
the excitation function above the 0.5 mb level. This 
two-phonon calculation is shown in Fig.|31as the dotted- 
dashed curve and it is evident that it does not reproduce 
the data well at low energies. We have also tried the 
recipe of a larger diffuseness (up to ai=5 fm) inside the 
Coulomb barrier, which was described in Ref. jj, but 
this did not result in any significant improvement of the 
fit to the data. 

A value of AR — 0.21 fm was needed to reproduce the 
present data in the cross section region of 0.1 to 100 /ib. 
The potential with this AR value produces a Coulomb 
barrier of 134 MeV, and a pocket inside the Coulomb 
barrier at 112.9 MeV, which is about 20 MeV higher than 
the ground state energy of the compound nucleus at 92.3 
MeV. The value of AR = 0.21 fm was also adopted in 
the 3ph calculations. The result is shown as the solid 
curve in Fig. \S\ which is seen to provide a better fit to 
the data in the 0.1 to 100 /ib cross section range. The 
dashed curve shows as a reference the results obtained in 
a one-dimensional calculation, i.e., without any couplings 
(with AR = 0.21 fm). It should be emphasized that the 
difference between the two sets of calculations, using the 
same AR = 0.21 fm, is rather modest. Hence, the 2ph 
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calculation with AR = 0.21 fm is not shown for clarity. 
These two calculations reproduce the experimental data 
in the highest energy region and in the 0.1 to 10 /zb range, 
but the cross sections in the region around E=130 MeV 
are overpredicted. These calculations also overpredict 
the cross sections at the lowest energies, which is the 
main topic of the present study. 

The coupled-channels calculations shown in Fig. |S1 ex- 
hibit essentially the same energy dependence (slope) as 
the one-dimensional calculation when the cross sections 
are small (< 10 /ib). The coupled- channels calculations 
are just shifted to lower energies relative to the one- 
dimensional calculation. In the present case the shift 
is about 7 MeV (for the same AR = 0.21 fm). This 
is a general feature of coupled-channels calculations and 
it is, therefore, very unlikely that any minor adjustment 
in the coupled-channels calculations would reproduce the 
steep falloff that the data exhibit at extreme sub-barrier 
energies. Thus, it appears that the fusion hindrance be- 
havior, which now has been observed for many systems, 
is also present in the new data for ^"^Ni-l-^^^Mo. This will 
be shown more convincingly in the next section, where 
other representations of the fusion cross section arc dis- 
cussed. 

From the experimental data and calculations alike, 
one can infer a rather broad, about 18 MeV wide, 
barrier distribution (not shown here). An elaborate 
coupling scheme is therefore needed in the coupled- 
channels calculations. At energies above the barrier, 
one-dimensional calculations result in higher cross sec- 
tions when compared to coupled-channels calculations. 
This behavior is caused by the long-range Coulomb- 
excitation/polarization of the low- lying quadrupole 
states, as pointed out in Ref. 0. As expected, the sup- 
pression of the coupled-channels fusion cross sections is 
rather strong in this region for the "soft" ^°°Mo nucleus. 
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FIG. 6: (a) Logarithmic derivative representation of the Ni 
-I- ^'"'Mo fusion excitation function, (b) S-factor represen- 
tation of the same data. All the circles were obtained from 
two successive data points and stars were derived from least 
square fits to three neighboring data points. Calculations are 
shown by curves. The extrapolation curve (dashed-dotted) 
shown in panel (b) was obtained from a straight line extrap- 
olation of the logarithmic derivative representation in panel 
(a). Included in (a) are also the data from Ref. X5] (solid 
triangles) and Ref. [T^ (open circles). See text for details. 



IV. LOGARITHMIC DERIVATIVES AND 
S-FACTORS 

The logarithmic derivative, L = d\n{aE)/dE, origi- 
nally introduced in Ref. Q , is shown in Fig. for the 
^■^Ni -I- -'^""Mo system. The solid circles were obtained di- 
rectly from two successive data points, whereas the stars 
were derived from least-squares fits to three neighboring 
data points. The lower limit of the logarithmic derivative 
was derived from the upper limit on the cross section at 
i;cm=120.2 MeV and the data point at i;cm=121.2 MeV 
(see Table 1). The present data are compared to those of 
Refs. 0,01, which are represented by open circles and 
solid triangles in Fig. respectively. Only two-point 
derivatives are shown for these data. The three-phonon 
coupled-channels and the one-dimensional barrier pene- 
tration calculations are shown as solid and dashed curves, 
respectively, while the thick solid, nearly horizontal line 
corresponds to a constant S'-factor expression derived in 
Ref. The present experimental results just reach 



the constant 5-factor line, implying that the experimen- 
tal data have reached a maximum value for the S'-factor. 
The energy, Es , representing the intersection between the 
experimental logarithmic derivatives and the constant 5- 
factor line corresponds to a value of £'^=120. 6 MeV. The 
dashed-dotted line in Fig.|H^ is an extrapolation obtained 
under the assumption that the logarithmic derivative is 
a straight line near the crossing point. This method was 
first introduced in Ref. ,10] in order to obtain the ex- 
trapolated values of Eg for the systems ^^Ni -|- ^°Ni and 
^*Ni -I- ^''Ni. The calculated logarithmic derivatives are 
seen to saturate around L=1.5 - 2 MeV~^ (or start to os- 
cillate) below i?cm= 124 MeV, whereas the corresponding 
experimental values continue to grow with decreasing en- 
ergies. This saturation behavior has already been noted 
in Refs. [Ill 113. 

The 5-factor representation for the *^'*Ni + ^°°Mo sys- 
tem is presented in Fig. ^jp. As the experimental loga- 
rithmic derivatives only just reach the constant S'-factor 
curve, the S-factor maximum is not fully developed. Ad- 
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TABLE IV: Parameters for Ni-induced fusion reactions on targets around the Z, N = 28, 40 and 50 shells. The parameters are: 
Z\Z2yfJi., where = A^A^ji^Ax + A2) is the reduced mass number of the coUiding system, the crossing point energy the 
lowest measured energy -Bmin, the reference energy B^'^^ obtained from Eq. the ratio EgjE^ , the minimum cross section 
measured, the fusion reaction Q-value, the Coulomb barrier (calculated with the Bass model fSOjl ') and the number of "valence 
nucleons", N^u, outside the nearest closed shells. For systems, which either exhibit a clear maximum in the 5- factor or whose 
logarithmic derivatives have not reached the constant S-factor curve yet, but can be extrapolated to obtain the crossing point, 
values Es are listed in column 3. For systems, where it is not possible to make good extrapolations, the lowest measured energy 
is quoted in column 4. For the Ni + Sn systems, where ^^^~^^*Sn targets have been used in the measurements, only data for 
the '^^''Sn target are included in the table for brevity. 
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ditional measurements would be required to clearly de- 
lineate the maximum in the S*- factor. The dashed-dotted 
curve corresponds to the straight line extrapolation of the 
logarithmic derivative in Fig. The two- and three- 
phonon calculations are shown in Fig. as dashed and 
solid curves, respectively. It is evident that coupled- 
channels calculations overpredict the fusion cross section 
at extreme sub-barrier energies. 



V. COMPARISON TO OTHER NI-INDUCED 
FUSION SYSTEMS 

The experimental fusion data involving Ni projec- 
tiles and compound nuclei in the A=100 - 200 region, 
are summarized in Table IV. For four systems, namely 
the 58Ni+58Ni, 64Ni+64Ni, 60]v^i_^89Y and the present 
^^Ni+^^^Mo reactions, the cross sections have been mea- 
sured to sufficiently low energies to determine the energy, 
Es of the maximum of the 5-factor representation. Pre- 
viously, we have obtained the Eg values for the ^*Ni+^°Ni 
and '^*Ni+^^Ni systems 9]. The values of E^ are listed 
in column three of Table IV. Two additiona l sy stems, 
namely ^SNi+^^Qe and ^'^Ni+^^'Ge from Ref.[ial, have 
been measured down to levels close to where the 5'-factor 
maximum occurs. The logarithmic derivative for these 
systems is shown in Fig. [T] The location of the S'-factor 



maximum was obtained by performing a small extrapo- 
lation (solid line) of the logarithmic derivative to where 
it crosses the constant 5*- factor line (dashed line). The 
resulting values are listed in Table IV. Because of the un- 
certainties inherent in these extrapolations, errors of 2% 
(~2 MeV) were assigned to the extrapolated Es data. Ta- 
ble IV also lists other fusion systems involving Ni beams 
for which the cross section was not measured down to a 
level that allows for determination of the S'-factor maxi- 
mum. For these systems only the lowest measured energy 
is listed. 

The values o f Eg are plotted a s a function of the pa- 
rameter ZiZ2^J A1A2/ {Ai + A2) in Fig.lHland compared 
to the empirical formula (solid curve), Eq. ^ obtained 
from a fit to all available fusion data involving stiff nuclei. 
Obviously, only two systems, ^'^Ni+'^'^Ni and 60Ni_|_89Y 
follow the systematics; all other systems fall below the 
curve. Previously, it has been pointed out [JOj that there 
is a rather compelling correlation between the stiffness of 
the interacting nuclei and the location of the S'-factor 
maximum, Es, relative to the empirical trend for the 
Ni+Ni systems. The addition to the systematics of the 
data points for *^'*Ni+i°°Mo and for the two Ni+Ge sys- 
tems, which all involve soft nuclei, appears to corrobo- 
rate this observation. Thus, for the ^^Ni+^°°Mo system, 
the value of El'^f predicted from Eq. [H is E^^^f =lM.h 
MeV whereas the measured value of £'^ = 120. 6 MeV is 
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only 90% of this value. The reason for this reduction 
is presumably that ^^Ni + ^'^"Mo should be viewed as 
an open-shell colliding system so that strong coupling ef- 
fects broaden the effective barrier distribution and push 
the energy where the hindrance behavior occurs down to 
even lower energies. 
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FIG. 7: Logarithmic derivative representations of the fusion 
data for the systems '^^Ni-h'^'^Ge (a) and ^''Ni+^'^Ge (b). The 
data are from Ref. [l^. 

The deviation of the measured or extrapolated values 
of Es from the expected E^^-^ value, seen in Fig. |S1 thus 
appears to depend on the stiffness of the fusing nuclei. 
A quantitative relation between the stiffness and the de- 
viation from E^'^-f is not yet known. As a first attempt, 
we associate the stiffness of a nucleus to its proximity 
to closed proton or neutron shells and define the num- 
ber of "valence nucleons", Np^, as the sum of particles 
and holes outside the nearest closed shells. Here, ^''Ni 
is considered to have four holes in the N=40 neutron 
shell rather than eight particles outside N=28. The val- 
ues of Nph are listed in Table IV and, in Fig.lHl the ratio 
Es/El^f is plotted as a function of this parameter. There 
is a general trend of decreasing values of Eg/E^'^^ with 
increasing values of Nph- We observe that for the other 
systems in Table IV, the data for E„iin (upper limits) 
are not in contradiction with Figs. |H1 and El ft should 
be noted that all of the systems shown in Figs. |S1 and |51 
and Table IV have rather large negative fusion Q— values. 
One may also compare the £^s-values to the height of the 
Coulomb barrier (obtained from the Bass prescription). 
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FIG. 8: Plot of vs. ZxZ^^ for Ni bombarding dif- 
ferent targets (see Table IV). Solid symbols correspond to 
systems for which the S'-factor maximum is well determined: 
5«Ni+58Ni, 64Ni+64Ni, «ONi+«»Y and ^^Ni+^ooMo. Open 
symbols are associated with the extrapolations for the systems 
58Ni+60Ni, 5»Ni-h'^*Ni, 58Ni-h^*Ge and ''"Ni-h^^Ge. The tri- 



angles represent Ni-|-Ni, diamonds Ni-fGe, a square 
and a circle ®'*Ni-|-^'"'Mo, respectively. 
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FIG. 9; Plot of Es/El^^ vs. Nph, where Npu is the total num- 
ber of "valence nucleons" outside closed shells in the entrance 
channel. Symbols are defined in Fig. |H1 



which is hstcd in Table IV. The ratio E./E^^f exhibit a 
stronger dependency on the value of Nph than that ob- 
tained for the ratio Eg/Vc indicating that the observed 
effect does not just depend on the change in the Coulomb 
barrier height with addition of neutrons to the interact- 
ing nuclei. 
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VI. DISCUSSION 

The above observations are all phenomcnological in na- 
ture as there is at present no satisfactory understanding 
of the fusion hindrance at extreme sub-barrier energies. 
Mariy authors have tried to explain this new phenomenon 
El IM El HI HI IM El- One suggestion is to use a 
large diffuseness of the ion-ion potential in the coupled- 
channels calculations. This recipe is sometimes used to 
explain high-precision fusion data It has been ar- 

gued that the failure to reproduce the steep fall-off is 
caused by the Hill- Wheeler approximation [^E 
Ref. i it was shown, however, that these limitations 
of the analysis are not responsible for the observed sub- 
barrier fusion hindrance phenomenon. Dasso and Pol- 
larolo tried to reduce the fusion cross sections by using 
a shallow well inside the barrier l3a . Giraud et al. in- 
vestigated the effect of "friction" [33. More discussions 
about this phenomenon can be found in Refs. E3) l4ll |. 
All of these suggestions may improve the fit of calcula- 
tions to the experimental data in some cases, but they 
do not provide a convincing explanation of the observed 
suppression in all systems. These studies are presently 
at the stage of raising questions and discovering weak 
points in the existing models. The data are most likely 
still insufficient to lead to the correct explanation. More 
precision sub-barrier fusion measurements are required 
to further explore which modifications in the theoreti- 
cal models are relevant for a correct description of the 
phenomenon. 

VII. SUMMARY AND CONCLUSION 

The phenomenon of sub-barrier fusion hindrance was 
first observed in systems involving stiff nuclei and a sim- 



ple expression J3, was derived for the energy at which 
the hindrance of fusion between such nuclei occurs. Fur- 
thermore, a study of Ni-|-Ni fusion involving different Ni 
isotopes has shown that the onset of the fusion hin- 
drance deviates from these systematics depending on the 
"stiffness" of the interacting nuclei. 

In the present work, we have measured the fusion ex- 
citation function for the system ^''Ni+^*'°Mo down to a 
cross section level of ~5 nb, i.e. about 12% below the 
Bass barrier. We observe that the fusion process is hin- 
dered at the lowest energies relative to expectations based 
on coupled-channels calculations. The present study of 
the ^^Ni-t-^^^Mo system shows that the fusion hindrance 
for this system, as well as two other soft systems, deviates 
strongly from the systematics, and it thus corroborates 
the earlier observation in Ni-|-Ni systems. It is further- 
more shown that this deviation depends monotonically 
on a parameter Nph, which is the sum of nucleons (holes) 
outside of closed shells of the fusing nuclei suggesting that 
this parameter is a good measure of the "stiffness" of the 
interacting nuclei. 

In conclusion, an interesting nuclear structure depen- 
dence of the fusion hindrance has been observed. The 
origin of this effect is still unknown. It occurs at rela- 
tively high excitation energies (for systems in Fig.|Hland 
El they are around Eex ~ 30 - 40 MeV) , where the natu- 
ral width of compound levels is larger than their spacing. 
A lack of available final states thus appears to be ruled 
out as an explanation. 
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